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Abstract 

In  this  paper,  we  exploit  the  notion  of  partial  sparsity  for  scene  reconstruction  associated  with  through-the-wall 
radar  imaging  of  stationary  targets  under  reduced  data  volume.  Partial  sparsity  implies  that  the  scene  being  imaged 
consists  of  a  sparse  part  and  a  dense  part,  with  the  support  of  the  latter  assumed  to  be  known.  For  the  problem  at 
hand,  sparsity  is  represented  by  a  few  stationary  indoor  targets,  whereas  the  high  scene  density  is  defined  by  exterior 
and  interior  walls.  Prior  knowledge  of  wall  positions  and  extent  may  be  available  either  through  building  blueprints  or 
from  prior  surveillance  operations.  The  contributions  of  the  exterior  and  interior  walls  are  removed  from  the  data 
through  the  use  of  projection  matrices,  which  are  determined  from  wall-  and  corner-specific  dictionaries.  The  projected 
data,  with  enhanced  sparsity,  is  then  processed  using  l}  norm  reconstruction  techniques.  Numerical  electromagnetic 
data  is  used  to  demonstrate  the  effectiveness  of  the  proposed  approach  for  imaging  stationary  indoor  scenes  using  a 
reduced  set  of  measurements. 
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1  Introduction 

The  ultimate  objective  of  achieving  actionable  intelligence 
in  an  efficient  and  reliable  manner  is  faced  with  a  host 
of  challenges  underlying  the  urban  sensing  and  through- 
the-wall  radar  imaging  (TWRI)  applications  [1-9].  First 
and  foremost  is  the  increasing  demand  on  radar  systems 
to  deliver  high-resolution  images  in  both  range  and 
crossrange,  which  requires  the  use  of  wideband  signals 
and  large  array  apertures,  respectively.  Second,  the 
backscatter  from  the  first  wall,  which  is  the  exterior 
wall  of  the  building  being  imaged,  is  much  stronger 
than  the  returns  from  the  interior  of  the  building.  This 
is  because  the  signal  undergoes  attenuation  in  the  wall 
materials  and  the  more  walls  the  signal  penetrates  to 
reach  the  indoor  targets,  the  weaker  are  the  returns. 
Therefore,  the  clutter  caused  by  the  wall  backscatter 
can  significantly  contaminate  the  radar  data  and  hinder 
the  main  intent  of  providing  enhanced  system  capabilities 
for  imaging  of  building  interiors  and  detection  and 
localization  of  stationary  indoor  targets. 

Recently,  compressive  sensing  (CS)  has  been  used  for 
efficient  data  acquisition  in  radar  systems  in  general 
[10-14]  and  in  urban  radar  systems  in  particular  [15-19]. 
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For  urban  radar  systems,  removal  of  clutter  and  stationary 
targets  via  change  detection  or  exploitation  of  sparsity  in 
the  Doppler  domain  readily  enables  the  application  of  CS 
techniques  for  moving  target  detection  inside  buildings 
[20-23].  However,  these  means  are  not  available  for  detec¬ 
tion  and  localization  of  stationary  targets  of  interest,  thereby 
requiring  significant  mitigation  of  wall  reflections. 

Several  approaches  have  been  proposed  for  dealing 
with  front  wall  returns  under  full  data  volume  without 
the  need  for  reference  or  background  data  [8,24-26].  In 
[8],  the  wall  parameters,  such  as  thickness  and  dielectric 
constant,  are  estimated  from  first  wave  arrivals  and  then 
used  to  model  and  subtract  the  front  wall  contributions 
from  the  received  data.  The  approach  in  [24]  for  wall 
clutter  mitigation  is  based  on  spatial  filtering.  It  utilizes 
the  strong  similarity  between  wall  electromagnetic  (EM) 
responses,  as  viewed  by  different  antennas  along  the  entire 
physical  or  synthesized  array  aperture.  A  spatial  filter  is 
applied  to  remove  the  dc  component  corresponding  to 
the  constant-type  radar  return  associated  with  the  front 
wall.  The  subspace  decomposition  method,  presented  in 
[25,26],  utilizes  not  only  the  approximately  identical  wall 
scattering  characteristics  across  the  array  elements  but 
also  the  higher  strength  of  the  wall  reflections  compared 
to  that  of  target  reflections.  When  singular  value  decom¬ 
position  (SVD)  is  applied  to  the  measured  data  matrix,  the 
wall  subspace  can  be  captured  by  the  singular  vectors 
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associated  with  the  dominant  singular  values.  As  a  result, 
the  wall  contributions  can  be  removed  by  projecting  the 
data  measurement  vector  at  each  antenna  on  the  wall  or¬ 
thogonal  subspace.  It  is  noted  that  as  the  round-trip  signal 
traveling  times  from  the  antennas  to  each  interior  wall, 
which  is  parallel  to  the  front  wall,  are  constant  across  the 
array  aperture,  both  spatial  filtering  and  the  subspace 
decomposition  methods  will  also  mitigate  returns  from 
interior  parallel  walls  as  long  as  they  are  not  shadowed 
by  other  contents  of  the  building  [27], 

Both  spatial  filtering  and  subspace  projection  approaches 
have  been  shown  to  be  equally  effective  for  synthetic  aper¬ 
ture  radar  (SAR)  imaging  under  reduced  data  volume, 
provided  that  the  same  reduced  set  of  frequencies  or  time 
samples  is  used  at  each  available  antenna  position  [28,29]. 
Requiring  the  same  frequency  observations  or  time 
samples  across  all  antennas  is  very  restrictive  and  may 
not  always  be  feasible.  For  example,  some  individual 
frequencies  or  frequency  subbands  may  be  unavailable 
due  to  competing  wireless  services  or  intentional  inter¬ 
ferences.  Further,  antenna  positions  may  signify  radar 
units  operating  independently,  each  with  a  separate 
frequency  band  to  avoid  cross-interference.  In  this 
paper,  we  propose  an  alternate  scheme  for  imaging  of 
stationary  indoor  scenes  which  overcomes  this  limitation 
of  wall  clutter  mitigation  techniques  under  reduced  data 
volume  by  exploiting  prior  knowledge  of  the  room  layout. 
This  information  may  be  available  either  through  building 
blueprints  or  from  prior  surveillance  operations  specific¬ 
ally  dedicated  to  determining  the  building  layout.  We  con¬ 
sider  the  scene  being  imaged  to  be  partially  sparse.  That 
is,  the  scene  consists  of  two  parts,  one  of  which  is  sparse 
and  contains  the  stationary  indoor  targets  of  interest  while 
the  other,  corresponding  to  the  exterior  and  interior  walls, 
is  dense  with  known  support  [30,31].  We  focus  on 
stepped-frequency  SAR  operation  and  assume  that  few 
frequency  observations  are  available,  which  could  be 
the  same  or  different  from  one  antenna  position  to  another, 
constituting  the  set  of  reduced  spatial  measurements.  We 
employ  projection  matrices  that  are  determined  from 
wall-  and  corner-specific  scattering  responses  to  remove 
the  exterior  and  interior  wall  contributions  from  the  mea¬ 
surements.  In  so  doing,  we  enable  the  application  of  con¬ 
ventional  sparse  reconstruction  schemes  to  obtain  an  image 
of  the  sparse  part  of  the  scene  containing  the  stationary  in¬ 
door  targets.  We  demonstrate  the  effectiveness  of  the  par¬ 
tial  sparsity-based  approach  for  reconstruction  of  stationary 
through-the-wall  scenes  using  numerical  electromagnetic 
data  of  a  single-story  building  for  both  cases  of  having  the 
same  reduced  set  of  frequencies  at  each  of  the  available  an¬ 
tenna  locations  and  also  when  different  frequency  measure¬ 
ments  are  employed  at  different  antenna  locations. 

It  is  noted  that,  as  an  alternative  to  the  proposed  ap¬ 
proach,  the  wall  and  corner  responses  can  be  modeled 


and  then  subtracted  from  the  received  data.  However, 
this  approach  has  two  issues.  First,  it  is  very  sensitive 
to  phase  errors.  Second,  it  may  not  always  work  since 
some  corners  and  parts  of  interior  walls  may  be  shadowed 
by  objects  in  the  interior  of  the  building.  Additionally,  it  is 
important  to  draw  a  distinction  between  the  proposed 
approach  and  the  subspace  projection  approach  for  wall 
removal  [25,26].  Although  both  approaches  involve  data 
projections,  they  exploit  fundamentally  different  charac¬ 
teristics  of  the  data  measurements  to  achieve  the  desired 
objective.  The  subspace  projection  approach  does  not  re¬ 
quire  knowledge  of  wall  positions  and  extent.  However,  it 
assumes  nominal  geometry  of  the  walls.  In  this  respect,  it 
relies  on  a  specific  radar  configuration,  whose  defining 
characteristics  are  the  constant  distance  of  each  antenna 
from  the  walls  and  normal  incidence  illuminations.  This 
ensures  approximately  constant  wall  contributions  in 
the  data  received  at  all  antennas.  The  partial  sparsity 
approach,  on  the  other  hand,  relies  on  accurate  knowledge 
of  the  building  geometry  to  create  wall-  and  corner- 
specific  dictionaries,  which  are  subsequently  used  for  miti¬ 
gating  the  contributions  of  exterior  and  interior  walls.  It 
does  not,  however,  demand  the  invariance  of  wall  returns 
across  the  partial  or  entire  radar  aperture.  The  impact  of 
these  fundamental  differences  on  the  performance  of  the 
two  approaches  is  highlighted  in  Section  5. 

The  remainder  of  this  paper  is  organized  as  follows. 
Section  2  presents  the  signal  model  under  the  assumption 
of  known  support  of  the  exterior  and  interior  walls.  The 
wall  contribution  removal  technique  and  scene  recon¬ 
struction  are  discussed  in  Section  3.  Section  4  evalu¬ 
ates  the  performance  of  the  proposed  partially  sparse 
through-the-wall  scene  reconstruction  approach  using 
numerical  EM  data  of  a  single-story  building.  Perform¬ 
ance  comparison  with  the  subspace  projection  ap¬ 
proach  is  also  provided.  Conclusions  are  drawn  in 
Section  5. 

2  Signal  model 

In  this  section,  we  develop  the  forward  scattering  model 
for  through-the-wall  radar  imaging.  The  model  is  based 
on  physical  optics  and  any  associated  nonlinearities  are 
ignored  [32]. 

Consider  a  monostatic  SAR  with  N  antenna  positions 
located  along  the  tf-axis  parallel  to  a  homogenous  front 
wall.  The  transmit  waveform  is  assumed  to  be  a  stepped- 
frequency  signal  of  M  frequencies,  which  are  equispaced 
over  the  desired  bandwidth  coM  ,  -  6J0, 

=  coq  +  mAca,  w  =  0,l,...,Af-l,  (1) 

where  cj0  is  the  lowest  frequency  in  the  desired  frequency 
band  and  A co  is  the  frequency  step  size.  The  scene  behind 
the  front  wall  is  assumed  to  be  composed  of  P  point 
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targets,  L- 1  interior  walls,  which  are  parallel  to  the  front 
wall  and  to  the  radar  scan  direction,  and  K  corners  cor¬ 
responding  to  the  junctions  of  two  walls  perpendicular 
to  each  other.  It  is  noted  that,  due  to  the  specular  nature 
of  the  wall  reflections,  a  SAR  system  located  parallel  to 
the  front  wall  will  only  be  able  to  receive  backscattered 
signals  from  interior  walls,  which  are  parallel  to  the 
front  wall.  The  contributions  of  walls  perpendicular  to 
the  front  wall  will  be  captured  primarily  through  the 
backscattered  signals  from  the  corners  [33,34], 

The  component  of  the  received  signal  corresponding 
to  the  mth  frequency  at  the  nth  antenna  position,  with 
phase  center  at  x,„  =  (xtn,  0),  due  to  the  P  point  targets  is 
given  by  [35,36] 

p- 1 

ztgt(m,  n)  =  exp(-jamTP'„) ,  (2) 

p= o 

where  ap  is  the  complex  amplitude  corresponding  to  the 
pth  target  return  and  Tpn  is  the  two-way  traveling  time 
between  the  nth  antenna  and  the  pth  target.  The  reflec¬ 
tions  from  the  L  walls  measured  at  the  nth  antenna  loca¬ 
tion  corresponding  to  the  mth  frequency  can  be  expressed 
as  [24] 

L- 1 

Zwaii(ra,  n)  =  exp  (-j(0mTwj) ,  (3) 

1=0 

where  aWti  is  the  complex  amplitude  associated  with  the 
/th  wall  and  twj  is  the  two-way  traveling  time  of  the  signal 
from  the  nth  antenna  to  the  /th  wall.  Note  that  since  the 
scan  direction  is  parallel  to  the  walls,  the  delay  twj  does 
not  depend  on  the  variable  n  and  is  a  function  only  of  the 
downrange  distance  between  the  /th  wall  and  the  antenna 
baseline.  Finally,  the  reflections  from  the  K  corners  mea¬ 
sured  at  the  nth  antenna  location  corresponding  to  the 
mth  frequency  can  be  expressed  as  [33,37] 

K- 1 

Zcorner  (m, «)  =  >]lT„o>smc  ( ( (omLk / c )  sin  (dk,„-8k))  exp 
k= 0 

(4) 

where  ak  is  the  complex  amplitude,  Lk  is  the  length,  0k 
is  the  orientation  angle  of  the  /<rth  corner,  ikn  is  the  two- 
way  propagation  delay  between  the  nth  antenna  and  the 
/(th  corner,  6kn  is  the  aspect  angle  associated  with  the 
/<th  corner  and  the  nth  antenna,  and  Tkt„  is  an  indicator 
function  which  assumes  a  unit  value  only  when  the  nth 
antenna  illuminates  the  concave  side  of  the  Ath  corner. 
We  note  that  each  of  the  complex  amplitudes  aPi  aWib  and 
d>  in  (2)  to  (4)  contains  contributions  from  free-space 
path  loss,  attenuation  due  to  propagation  through  the  wall 
(s),  and  the  reflectivity  of  the  corresponding  scatterer.  The 


nth  received  signal  corresponding  to  the  mth  frequency  is, 
thus,  given  by 

z(m,  n)  =  Zwaii(m,  n)  +  zcorner(m,  n)  +  ztgt(m,  n)  (5) 

Assume  that  the  scene  being  imaged  is  divided  into  a 
finite  number  of  grid-points,  say  Q,  in  crossrange  and 
downrange.  Let  z„  represent  the  received  signal  vector 
corresponding  to  the  M  frequencies  and  the  nth  antenna 
location,  and  s  be  the  concatenated  Q  x  1  scene  reflectivity 
vector  corresponding  to  the  spatial  sampling  grid.  Under 
the  assumption  that  the  building  layout  is  known  a  priori,  s 

can  be  expressed  as  =  [sf  sj] '  ,  where  S]  eCQl  is  the  dense 
part  whose  support  is  known  and  s2sC^2,  Q2  =  Q-Qi,  is 
the  sparse  part.  Note  that  corresponds  to  the  walls 
that  are  parallel  to  the  antenna  baseline.  Further,  since 
the  wall  junctions  lie  along  the  parallel  walls,  the  corner 
locations  would  correspond  to  the  support  of  a  subset 
of  Si,  say  SieCR,R  <  Q, .  Then,  using  (2)  to  (5),  we  obtain 
the  matrix -vector  form 

z„  —  AnS\  T  BnS\  ~F  Cns2 ,  (6) 

where  An,  B„,  and  C„  are  the  dictionary  matrices  corre¬ 
sponding  to  the  wall,  corner  reflector,  and  point  target, 
respectively.  The  matrix  C„  is  of  dimension  Mx  Q2  with 
its  (m,  qo)th  element  given  by 

\Cn\m,q2  =  exp(-/«mr?2j„)  (7) 

where  is  the  two-way  traveling  time  between  the  nth 
antenna  and  the  q2 th  grid-point  of  the  sparse  part.  The 
wall  dictionary  An  is  an  M  x  Qx  matrix,  whose  (m,  qk) th 
element  takes  the  form  [38] 

\An\m,qi  =  exp  (-;wm2y?i/c)3?lK  (8) 

In  (8),  y  represents  the  downrange  coordinate  of  the 
<7ith  grid-point  in  the  dense  part,  and  3?|  M  is  an  indicator 
function,  which  assumes  a  unit  value  only  when  the  qk th 
grid-point  lies  in  front  of  the  nth  antenna,  as  illustrated 
in  Figure  1.  That  is,  if  xqi  represents  the  crossrange  co¬ 
ordinate  of  the  r/ith  dense  grid-point  and  Sx  represents 
the  crossrange  sampling  step,  then  3?|  „  =  1  provided 
that  xqi~Y  - xtn  ^  xql  +  2  •  corner  dictionary  Bn  is 
anMxi?  matrix  whose  (m,  r)th  element  is  given  by 

[Bn]mr  =  exp(-/wmrr,„)r,.i„sinc^6JyTj  sin (0r,H  - 0,-) J 

(9) 

Equation  6  considers  the  contribution  of  only  one 
antenna  location.  Stacking  the  measurement  vectors 
corresponding  to  all  N  antennas  to  form  a  tall  vector, 
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qqt/i  pixel 


U 

I  %tn 

Figure  1  Illustration  of  the  indicator  function  for  the  wall 
returns.  It  depicts  that  the  indicator  function  will  assume  unit 


values  only  for  the  gray  pixels  when  the  antenna  is  at  the 
position  shown. 

r  t  t 

Z  =  1*1  z2 

zT  Y 

zN-l\  5 

(10) 

we  obtain  the  linear  system  of  equations 

z  =  Asi  +  Bsi 

+  CS2 , 

(11) 

where 

E-l  rH 

Sh  o 

II 

(12) 

B  =  [Btq  Bf  • 

c  =  [c0r  cf 

-Ctn_x]T 

The  vector  z  contains  the  full  dataset  corresponding  to 
the  N  antenna  locations  and  the  M  frequencies.  For  the 
case  of  reduced  data  volume,  consider  f),  which  is  a  / 


Table  1  Complex  dielectric  constant  for  materials  used  in 
the  simulation  code 


Material 

£' 

e" 

Brick 

3.8 

0.24 

Concrete 

6.8 

1.2 

Glass 

6.4 

0 

Wood 

2.5 

0.05 

Sheetrock 

2.0 

0 

Foam  cushion  and  fabric 

1.4 

0 

Ground 

10 

0.6 

Human 

50 

12 

(<</WA/)-dimensional  vector  consisting  of  elements  ran¬ 
domly  chosen  from  z  as  follows: 

^  =  Oz  =  OAsi  T  OBsr  T  OCS2  (13) 

In  (13),  ®  is  a  /  x  MN  measurement  matrix  of  the  form 

®  =  kron(^,//1)  ■  diag{#>0,  </)l, ..., 

/  =  /1/2, 

(14) 

where  ‘kron’  denotes  the  Kronecker  product,  //,  is  a 
/,  x  /,  identity  matrix,  i//  is  a  /2  x  N  measurement  matrix 
constructed  by  randomly  selecting  /2  rows  of  an  N  x  N 
identity  matrix,  and  (pm  n  =  0,  1,...,  N  -  1,  is  a  Jx  x  M 
measurement  matrix  constructed  by  randomly  selecting 
Ji  rows  of  an  M  x  M  identity  matrix.  We  note  that  xj/ 
determines  the  reduced  antenna  locations,  whereas  (f>n 
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determines  the  reduced  set  of  frequencies  corresponding 
to  the  «th  antenna  location. 

3  Wall  contribution  removal  and  scene 
reconstruction 

Given  the  reduced  measurement  vector  £;  and  knowledge 
of  the  support  of  the  walls  and  corners,  the  goal  is  to  re¬ 
construct  the  sparse  part  of  the  image  where  the  targets 
of  interest  are  located.  Towards  this  goal,  we  first  need 
to  remove  the  contributions  of  the  dense  part  of  the 
scene  from  Let  PA  be  the  matrix  of  the  orthogonal 
projection  from  CQ  onto  the  orthogonal  complement  of 


the  range  space  of  the  matrix  ®A  If  A>A  is  a  full  rank 
matrix,  then  PA  can  be  expressed  as 

Pa  =  I,-{OA)(OA)f  (15) 

where  Ij  is  a  /  x  /  identity  matrix  and  (®A)+  denotes  the 
pseudoinverse  of  (®^4).  On  the  other  hand,  if  OA  has  a 
reduced  rank,  then  we  have  to  resort  to  the  SVD  of  <S>A 
to  obtain  the  matrix  PA  as 

PA  =  UaUha  (16) 

where  UA  is  the  matrix  consisting  of  the  left  singular 
vectors  corresponding  to  the  zero  singular  values  and 


(b) 

Figure  3  Backprojection  results  using  raw  data  under  the  full  data  volume,  (a)  Full  image,  (b)  Image  with  the  support  region  of  the  exterior 
and  interior  walls  masked  out. 
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the  superscript  ‘Pf  denotes  the  Hermitian  operation.  Ap¬ 
plying  the  projection  matrix  PA  to  the  observation  vector  £j, 
we  obtain 

£,a=PaK  =  PA{<t>ASl  +  +  0 >Cs2)  (17) 

=  PA(t>Bs}  +  PA<l>C$2 

Next,  consider  the  projection  matrix  PB  given  by 

p  j  I,-(PA<]>B)(PA<t>B?  \{PA<t>B  has  a  full  rank 

B  \  UbUb  otherwise 

(18) 

where  11  b  is  the  matrix  consisting  of  the  left  singular 
vectors  corresponding  to  the  zero  singular  values  of  the 


matrix  PAOB.  Application  of  PB  to  the  measurement 
vector  S/t  leads  to 

Sba=Pb£,a  =  PbJPa^Bs]  +  PA®Cs2)  (19) 

=  PbPa(PCs2 

Thus,  after  the  sequential  application  of  the  two  projec¬ 
tion  matrices,  the  measurement  vector  ^BA  contains  con¬ 
tributions  from  only  the  sparse  image  part,  s2,  which  can 
then  be  recovered  by  solving  the  problem 

s2=  arg  ||s2||iSubjectto?M  ~PBPAOCs2  (20) 

The  problem  in  (20)  belongs  to  the  classical  setting  of  CS 
and,  thus,  can  be  solved  using  convex  relaxation,  greedy 


(b) 

Figure  4  Backprojection  results  after  application  of  the  subspace  projection  technique  under  the  full  data  volume,  (a)  Full  image, 
(b)  Image  with  the  support  region  of  the  exterior  and  interior  walls  masked  out. 
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pursuit,  or  combinatorial  algorithms  [39-43].  In  this  work, 
we  choose  orthogonal  matching  pursuit  (OMP),  which  is 
an  iterative  greedy  algorithm  [44]. 

It  is  noted  that  the  dimensionality  of  the  orthogonal 
complements  of  the  range  spaces  of  the  matrices  ®A 
and  ®B  is  at  least  /-  Qi  and  J  -  R,  respectively.  Further, 
if  access  to  full  data  volume  is  available,  the  proposed 
wall  removal  procedure  can  also  be  applied  to  the  full  data 
vector  z  with  appropriate  projection  matrices  determined 
from  the  dictionary  matrices  A  and  B  instead  of  ®A  and 
OB,  respectively.  The  wall-free  data  can  then  be  processed 
using  conventional  image  formation  techniques  [45] . 

4  Simulation  results 

In  this  section,  we  present  scene  reconstruction  results 
for  the  partial  sparsity  technique  using  numerical  EM  data 
and  provide  performance  comparison  with  the  subspace 
projection  approach  for  both  full  and  reduced  data  vol¬ 
umes.  For  the  reduced  data  volume,  we  consider  both 
cases  of  having  different  frequency  measurements  at  dif¬ 
ferent  available  antenna  locations  and  also  when  the  same 
reduced  set  of  frequencies  is  employed  at  each  of  the 
available  antenna  locations.  Note  that  for  the  subspace 
projection-based  wall  mitigation  CS  approach  proposed  in 
[28],  the  former  casts  a  more  challenging  problem  than 
the  latter,  as  it  is  not  amenable  to  wall  removal  using  dir¬ 
ect  implementation  of  the  subspace  projection  technique. 
Instead,  the  range  profile  at  each  employed  antenna  loca¬ 
tion  first  needs  to  be  reconstructed  through  /,  norm 
minimization  using  the  reduced  frequency  set  [28].  Then, 
the  Fourier  transform  of  each  reconstructed  range  profile 
is  taken  to  recover  the  full  frequency  data  measurements 
at  each  antenna  location.  Direct  application  of  the  sub¬ 
space  projection  technique  can  then  proceed,  followed  by 
the  scene  reconstruction. 


Table  2  TCR  -  backprojection  under  full  data  volume 


Method 

Figure  number 

TCR  (dB) 

Raw  data 

3b 

-4.05 

Subspace  projection 

4b 

-3.97 

Proposed 

5 

-1.45 

4.1  Electromagnetic  modeling 

The  simulation  is  based  on  XpatclT,  developed  by 
SAIC/DEMACO  (Champaign,  IL,  USA),  which  is  a  com¬ 
putational  EM  code  implementing  an  approximate  ray  tra¬ 
cing/physical  optics  computational  approach.  We  created 
the  computer  model  of  a  single-story  building,  with  over¬ 
all  dimensions  of  7  m  x  10  m  x  2.2  m,  containing  four 
humans  (labeled  1  through  4)  and  several  furniture  items, 
as  shown  in  Figure  2.  The  origin  of  the  coordinate  system 
was  chosen  to  be  in  the  center  of  the  building,  with  the 
x-axis  and  the  y-axis  oriented  as  shown  in  Figure  2b. 
The  exterior  walls  were  made  of  0.2-m-thick  bricks  and 
had  glass  windows  and  a  wooden  door.  The  interior  walls 
were  made  of  5-cm-thick  Sheetrock  and  had  a  wooden 
door.  The  ceiling/roof  is  flat,  made  of  a  7. 5-cm-thick 
concrete  slab.  The  entire  building  is  placed  on  top  of  a 
dielectric  ground  plane.  The  furniture  items,  namely,  a 
bed,  a  couch,  a  bookshelf,  a  dresser,  and  a  table  with 
four  chairs,  were  made  of  wood,  while  the  mattress  and 
cushions  were  made  of  a  generic  foam/fabric  material. 
Humans  1  through  4  were  positioned  at  various  locations 
in  the  interior  of  the  building  with  45°,  0°,  -20°,  and  10° 
azimuthal  orientation  angles,  respectively.  Note  that  an 
orientation  angle  of  0°  corresponds  to  the  human  facing 
along  the  positive  x  direction  and  the  positive  angles  cor¬ 
respond  to  a  counterclockwise  rotation  in  the  horizontal 
plane.  Human  3,  positioned  inside  the  interior  room,  was 
carrying  an  AK-47  rifle.  The  human  model  was  made  of  a 
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uniform  dielectric  material  with  properties  close  to  those 
of  the  skin  and  is  described  in  [46] .  The  human  body  radar 
cross  section  (RCS)  depends  on  the  aspect  angle  but  is 
generally  bounded  between  -10  and  0  dBsm.  Interestingly, 
the  average  human  body  RCS  is  fairly  constant  over  the 
frequency  range  considered  in  this  paper.  More  detailed 
results  on  the  human  body  radar  signature  can  be  found 
in  [46] .  The  AK-47  model  is  made  of  metal  and  wood  and 
is  described  in  [47,48].  The  dielectric  properties  of  the 
various  materials  employed  are  listed  in  Table  1. 

A  6-m-long  synthetic  aperture  line  array,  with  an 
inter-element  spacing  of  2.54  cm  and  located  parallel  to 
the  front  of  the  building  at  a  standoff  distance  of  4  m, 
was  used  for  data  collection.  Monostatic  operation  was 
assumed.  The  antenna  had  a  3-dB  beamwidth  of  60°  in 
both  elevation  and  azimuth  and  was  positioned  0.5  m 
above  the  ground  plane.  The  antenna  boresight  was  aimed 
perpendicular  to  the  exterior  wall.  A  stepped-frequency 
signal  covering  the  0.7-  to  2-GHz  frequency  band  with  a 
step  size  of  8.79  MHz  was  employed.  Thus,  at  each  of  the 
239  scan  positions,  the  radar  collected  148  frequency  mea¬ 
surements  over  the  1.3-GHz  bandwidth. 

4.2  Image  reconstruction  under  full  data  volume 

The  region  to  be  imaged  was  chosen  to  be  9  m  x  12  m 
centered  at  the  origin  and  divided  into  121  x  161  pixels, 
respectively.  Figure  3a  shows  the  image  obtained  with 
backprojection  using  the  full  raw  dataset.  In  this  figure 
and  all  subsequent  figures  in  this  paper,  we  plot  the 
image  intensity  with  the  maximum  intensity  value  in 
each  image  normalized  to  0  dB.  The  Hanning  window 
was  applied  to  the  data  along  the  frequency  dimension 
in  order  to  reduce  the  range  sidelobes  in  the  image.  The 
humans  in  the  image  are  indicated  by  red  circles.  We 
can  clearly  see  the  front  wall,  some  of  the  corners,  and 


Figure  7  An  illustration  of  the  random  subsampling  pattern 
employed  for  data  reduction.  Different  sets  of  reduced 
frequencies  are  employed  at  each  antenna  location. 


humans  1  and  2.  Human  3  in  the  interior  room  is  barely 
visible  due  to  the  additional  EM  loss  as  the  transmitted 
signal  has  to  penetrate  through  both  the  exterior  and 
interior  walls.  Likewise,  it  is  a  challenge  to  detect  human 
4,  who  is  the  farthest  away  from  the  front  wall.  Figure  3b 
shows  the  backprojected  image  after  masking  out  the 
dense  regions  with  known  support.  Although  all  the  tar¬ 
gets  are  visible  in  the  masked  image,  the  image  is  cluttered 
due  to  the  presence  of  the  wall  and  corner  sidelobes. 

Next,  we  reconstructed  the  scene  using  the  subspace 
decomposition-based  wall  mitigation  approach.  The  first 
two  dominant  singular  vectors  of  the  frequency  vs.  antenna 
raw  data  matrix  were  used  to  reconstruct  the  wall  subspace. 
The  wall  subspace  dimension  was  selected  using  the 
method  reported  in  [49].  Finally,  backprojection  was 
performed  on  the  wall  clutter  mitigated  data,  and  the 
corresponding  image  is  shown  in  Figure  4a.  We  observe 
that  although  the  stationary  targets  are  more  visible  and 
the  front  and  interior  wall  reflections  are  successfully 
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removed,  the  corners  indicating  the  presence  of  doors  and 
windows  are  still  present.  So  is  most  of  the  back  wall  due 
to  shadowing  effects.  The  approach  also  removed  the 
reflections  from  the  edge  of  the  couch,  and  only  the 
couch  corners  survive.  More  importantly,  the  presence 
of  discontinuities  in  the  front  wall  (windows  and  door) 
causes  the  subspace  decomposition-based  approach  to 
introduce  artifacts  in  the  image,  indicated  by  the  red 
rectangles.  Such  artifacts  in  the  interior  of  the  building 
are  more  visible  in  Figure  4b,  which  shows  the  image  after 
masking  out  the  dense  regions  with  known  support.  These 
artifacts  are  attributed  to  the  fact  that  the  subspace 


projection  scheme  assumes  the  wall  response  to  be  the 
same  from  one  antenna  to  the  other,  which  is  violated  by 
the  presence  of  windows  and  doors. 

Finally,  Figure  5  presents  the  backprojected  image  ob¬ 
tained  using  the  proposed  approach.  The  dense  part  of  the 
scene,  corresponding  to  the  building  layout  (exterior  and 
interior  walls  parallel  to  the  array  and  corners),  consisted 
of  7,196  pixels,  while  the  sparse  part  of  the  scene  consisted 
of  the  remaining  12,285  pixels.  Compared  to  Figures  3b 
and  4b,  the  image  in  Figure  5  is  the  least  cluttered  since 
the  wall  sidelobes,  in  particular  near  the  back  wall,  are  ab¬ 
sent.  All  of  the  humans  and  the  furniture  items  are  clearly 


(b) 

Figure  8  Scene  reconstruction  corresponding  to  the  subsampling  in  Figure  7.  (a)  Partial  sparsity-based  approach,  (b)  Subspace  projection 
wall  mitigation-based  CS  approach. 
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Table  3  TCR  -  sparse  reconstruction  when  different  sets 
of  reduced  frequencies  are  used  from  different  antenna 
locations 


Method 

Figure  number 

TCR  (dB) 

Subspace  projection 

8a 

2.75 

Proposed  (sparsity  =  10) 

8b 

14.96 

Proposed  (sparsity  =  25) 

11a 

10.12 

Proposed  (sparsity  =  45) 

11b 

6.86 

Proposed  with  cross-validation 

12 

24.78 

visible  in  the  image.  We,  therefore,  conclude  that  the 
proposed  approach  provides  superior  performance  com¬ 
pared  to  the  subspace  decomposition-based  wall  mitiga¬ 
tion  approach  under  the  full  data  volume. 

In  addition  to  the  visual  inspections,  we  also  assess  the 
performance  of  the  various  methods  in  terms  of  the 
target-to-clutter  ratio  (TCR),  which  is  defined  as  the  ra¬ 
tio  between  the  average  pixel  power  It  in  the  target  re¬ 
gion  to  the  average  pixel  power  Ic  in  the  clutter  region 
of  the  reconstructed  image  s2 

TCR=101og10(JY  (21) 


where  Rt  is  the  target  region,  R,  is  the  clutter  area,  Nt  is 
the  number  of  pixels  in  the  target  area,  and  Nc  is  the 
number  of  pixels  in  the  clutter  region.  The  target  area  is 
composed  of  regions  containing  the  four  humans,  and 
the  remaining  pixels  of  s2  constitute  the  clutter  region. 


Note  that  we  consider  furniture  reflections  as  unwanted 
returns,  and  accordingly,  they  are  treated  as  clutter.  Table  2 
shows  the  TCR  values  for  the  reconstructed  images  of 
Figures  3b,  4b,  and  5.  As  expected,  the  TCR  is  improved 
using  the  proposed  scheme  over  the  subspace  projection- 
based  wall  clutter  mitigation  scheme. 

4.3  Image  reconstruction  under  different  sets  of  reduced 
frequencies  at  each  available  antenna  location 

Conventional  image  formation  techniques,  such  as  back- 
projection,  compromise  the  image  quality  when  a  reduced 
number  of  measurements  is  considered,  thereby  impeding 
the  detection  of  targets  behind  the  wall  in  the  image 
domain.  This  is  illustrated  in  Figure  6,  wherein  we  used 
118  randomly  selected  frequencies  (79.7%  of  148)  and  79 
randomly  chosen  antenna  locations  (30%  of  239)  for 
backprojection,  which  collectively  represent  26%  of  the 
total  data  volume.  The  corresponding  space-frequency 
sampling  pattern  is  shown  in  Figure  7,  where  the  verti¬ 
cal  axis  represents  the  antenna  location  and  the  hori¬ 
zontal  axis  represents  the  frequency.  The  filled  boxes 
represent  the  data  samples  constituting  the  reduced  set 
of  measurements. 

Next,  we  reconstructed  the  scene  using  the  partial 
sparsity  approach  with  26%  data  volume.  The  number 
of  OMP  iterations,  usually  associated  with  the  sparsity 
level  of  the  scene,  was  set  to  10.  In  this  case,  and  for  all 
subsequent  sparse  imaging  results,  each  imaged  pixel  is 
the  result  of  averaging  200  runs,  with  a  different  random 
selection  for  each  run.  The  partial  sparsity-based  re¬ 
construction  of  the  sparse  part  of  the  scene  is  shown  in 
Figure  8a.  We  observe  from  Figure  8a  that  the  partial 
sparsity-based  scheme  was  able  to  detect  and  localize 
humans  1  through  3  successfully,  while  it  missed  human 
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Figure  9  An  illustration  of  the  random  subsampling  pattern  employed  for  data  reduction.  Same  set  of  reduced  frequencies  is  employed 

from  different  antenna  locations. 
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4.  In  addition,  some  clutter  (arising  from  the  left  chair  and 
table)  and  background  noise  is  also  visible  in  the  recon¬ 
structed  image. 

For  comparison,  we  also  performed  scene  reconstruc¬ 
tion  using  the  subspace  projection-based  wall  mitigation 
CS  approach  of  [28]  with  26%  data  volume.  Full  frequency 
data  measurements  were  first  recovered  from  the  lx  norm- 
reconstructed  range  profiles  at  each  considered  antenna 
location.  The  number  of  OMP  iterations  was  set  to  100 
for  each  range  profile  reconstruction.  This  is  because  the 
presence  of  the  wall  return  and  clutter  renders  the  range 
profile  quite  dense.  The  subspace  projection  approach 
was  then  applied  wherein  the  first  two  dominant  singular 


Table  4  TCR  -  sparse  reconstruction  when  the  same  set  of 
reduced  frequencies  is  used  from  different  antennas 


Method 

Figure  number 

TCR  (dB) 

Subspace  projection 

10a 

12.25 

Proposed 

10b 

13.32 

vectors  of  the  148  x  79  data  matrix  were  used  to  recon¬ 
struct  the  wall  subspace.  Finally,  standard  lx  norm  image 
reconstruction  was  performed  on  the  wall  clutter  mitigated 
full  frequency  recovered  data  to  form  an  image  of  the 
sparse  part  of  the  scene,  shown  in  Figure  8b.  Similar  to  the 
partial  sparsity  approach,  the  number  of  OMP  iterations  in 


(b) 

Figure  10  Scene  reconstruction  corresponding  to  the  subsampling  in  Figure  9.  (a)  Partial  sparsity-based  approach  (b)  Subspace  projection 
wall  mitigation-based  CS  approach. 
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this  case  was  chosen  to  be  10.  We  observe  from  Figure  5b 
that  human  1  was  detected,  human  2  was  barely  detected, 
while  humans  3  and  4  were  both  missing  from  the  recon¬ 
struction.  Moreover,  significantly  more  clutter  and  noise 
was  reconstructed  compared  to  Figure  8a.  We,  therefore, 
conclude  that  the  partial  sparsity-based  approach  com¬ 
pared  to  the  subspace  projection-based  wall  mitigation  CS 
approach  provides  superior  performance  for  the  same  re¬ 
duced  data  volume  when  different  sets  of  frequencies  are 
employed  at  the  available  antennas.  This  is  also  confirmed 
by  the  corresponding  TCR  values  provided  in  Table  3  (first 
two  rows). 


4.4  Image  reconstruction  under  the  same  set  of  reduced 
frequencies  at  each  available  antenna  location 

We  now  proceed  with  image  reconstruction  when  the 
same  reduced  set  of  frequencies  is  employed  at  each  of 
the  available  antenna  locations.  The  corresponding  space- 
frequency  sampling  pattern  is  shown  in  Figure  9.  We  use 
the  same  set  of  118  randomly  selected  frequencies  (79.7% 
of  148)  at  each  of  the  79  randomly  chosen  antenna  loca¬ 
tions  (30%  of  239).  Figure  10a  presents  the  result  of  the 
partial  sparsity-based  approach  with  the  number  of  OMP 
iterations  set  to  10.  We  observe  from  Figure  10a  that 
humans  1  through  3  were  successfully  localized,  while 


(b) 

Figure  1 1  Scene  reconstruction  using  partial  sparsity-based  approach  with  (a)  25  and  (b)  45  OMP  iterations. 

s _ I _ / 
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human  4  was  missed.  In  addition,  some  clutter  arising 
from  the  furniture  and  noise  were  also  reconstructed. 

We  next  applied  the  subspace  projection-based  wall 
mitigation  approach  directly  to  the  reduced  dimension 
data  matrix,  118  x  79  instead  of  148  x  239  [28].  The  wall 
suppressed  data  was  then  used  to  obtain  the  /,  norm- 
reconstructed  image  of  the  sparse  part  of  the  scene,  shown 
in  Figure  10b,  obtained  using  OMP  with  10  iterations.  We 
observe  that  the  subspace  projection  scheme  was  able  to 
detect  and  localize  humans  1,  2,  and  4  successfully,  while 
human  3  was  barely  detected.  Residual  wall  clutter  and 
some  of  the  furniture  returns  are  also  visible  in  the 
reconstructed  image.  We,  therefore,  conclude  that  the 
partial  sparsity-  and  the  subspace  decomposition-based 
wall  mitigation  CS  approaches  provide  comparable  per¬ 
formance  for  the  same  reduced  data  volume  when  the 
same  set  of  frequencies  is  employed  at  the  available  anten¬ 
nas.  This  is  validated  by  Table  4  which  shows  that  the  two 
schemes  have  comparable  TCRs. 

4.5  A  note  on  the  number  of  OMP  iterations 

OMP,  like  other  greedy  iterative  algorithms,  requires  the 
specification  of  the  scene  sparsity  for  exact  reconstruction 
[44].  In  most  practical  situations,  including  through-the- 
wall  imaging,  this  information  is  not  available  a  priori. 
Therefore,  the  stopping  criterion  based  on  the  fixed  num¬ 
ber  of  iterations,  which  is  tied  to  the  scene  sparsity,  is 
heuristic.  Figure  lla,b  shows  the  reconstruction  result  for 
the  scene  in  Figure  2  using  26%  of  the  data  volume  with 
the  number  of  OMP  iterations  chosen  to  be  25  and  45, 
respectively.  The  space-frequency  sampling  pattern  of 
Figure  7  was  employed.  In  both  cases,  humans  1  through 
3  are  clearly  visible.  However,  a  higher  amount  of  clutter 
and  background  noise  is  reconstructed  with  increasing 


number  of  iterations.  The  corresponding  TCR  values  are 
provided  in  Table  3  (rows  3  and  4). 

Use  of  cross-validation  has  been  proposed  to  prevent 
early/late  termination  of  greedy  reconstruction  algorithms 
[50].  Cross-validation  is  a  statistical  technique  that 
separates  a  dataset  into  a  training/estimation  set  and  a 
test/cross-validation  set.  The  test  set  is  used  to  prevent 
underfitting/overfitting  on  the  training  set.  The  cross¬ 
validation-based  OMP  reconstruction  result  using  26% 
of  the  data  volume  is  depicted  in  Figure  12,  with  one 
fifth  of  the  measurements  used  for  cross-validation. 
We  observe  that  the  cross-validation-based  approach 
fails  to  solve  the  problem  and  only  humans  1  and  2  are 
visible  in  the  reconstructed  image.  This  is  because  the 
signal  strength  from  humans  3  and  4  is  either  comparable 
to  or  weaker  than  those  from  sources  of  clutter,  and  the 
cross-validation-based  approach  regards  humans  3  and  4 
as  part  of  the  background  noise  level  [51].  Although  only 
two  of  the  four  humans  are  detected,  the  corresponding 
TCR  value  is  quite  high,  as  shown  in  Table  3  (last  row). 
This  is  because  very  little  clutter  is  reconstructed.  Various 
adaptive  approaches  have  recently  been  proposed  to  coun¬ 
ter  the  problem  of  low  signal-to-noise-and-clutter  ratio 
[51,52].  The  offering  of  these  schemes  to  the  problem  at 
hand  remains  to  be  explored. 

5  Conclusions 

In  this  paper,  we  applied  partial  sparsity  to  scene  recon¬ 
struction  associated  with  through-the-wall  radar  imaging 
of  stationary  targets.  Partially  sparse  recovery  deals  with 
the  case  when  it  is  known  a  priori  that  part  of  the  scene 
being  imaged  is  dense  while  the  rest  is  sparse.  For  the 
underlying  problem,  the  dense  part  of  the  scene  corre¬ 
sponds  to  the  building  layout  and  the  support  of  the 
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Figure  12  Partially  sparse  scene  reconstruction  using  cross-validation-based  OMP. 
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corresponding  part  of  the  image  is  assumed  to  be  known 
beforehand.  This  knowledge  may  be  available  either 
through  building  blueprints  or  from  prior  surveillance 
operations.  Using  numerical  EM  data  of  a  single-story 
building,  we  demonstrated  the  effectiveness  of  the  par¬ 
tially  sparse  reconstruction  in  detecting  and  locating 
stationary  targets  in  through-the-wall  scenes  while 
achieving  a  sizable  reduction  in  the  data  volume. 

Abbreviations 

CS:  compressive  sensing;  EM:  electromagnetic;  OMP:  orthogonal  matching 
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